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I. Introduction 

This report contains user instructions, a listing and 
documentation for a microcomputer BASIC program that can be used 
to compute an estimate of the probability that an electrometer 
based detection system will detect a submarine electric dipole 
field during an encounter. 

The program generates detection probabilities based on two 
encounter models. In the first encounter model, the detection 
system uses a crosscorrelation detector. In the second encounter 
model, the detection system uses a square law detector. Relative 
to operationally realizable values, probabilities based on the 
first model represent upper bounds and those based on the second 
model represent lower bounds. For both encounter models, the 
signal is proportional to the square of the component of the 
electric field that is measured by a single axis electrometer in 
an infinite medium of constant conductivity and the electric 
field noise does not change with changes in the position of the 
electrometer. Also, for both encounter models, the encounter is 
a straight line encounter with constant vertical separation. 

The encounter models can be interpreted as models of a 
electric anomaly detection system that is moving with constant 
course, speed and altitude in an encounter with a submarine 
moving with constant course, speed and depth. Or, they can be 
interpreted as models of a stationary electric anomaly detection 
system in an encounter with a submarine moving with constant 
course, speed and depth. 
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The program parameters include: electrometer axis 

direction, electrometer axis depression angle, encounter medium 
conductivity, submarine electric moments, submarine course, speed 
and altitude, electrometer course, speed and depth, encounter 
lateral range (the horizontal range at the closest point of 
approach in a straight line encounter) and false alarm rate. (In 
the program, a false alarm is the event that the detection system 
classifies noise as a dipole signal.) 
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II. Program User Instructions 

As listed in Appendix 10, the program can be run under a 
BASIC language that is compatible with IBM PC BASIC. If the 
listing is used to enter the program through the keyboard, then 
the program should be saved with the name EAD.BAS or the value of 
N$ on line 40 should be changed to the file name under which it 
is saved. 

The program contains user instructions in the form of query 
and parameter limitation messages. As an example of the former, 
after starting the program under BASIC, the following message 
should appear: 

Electric Anomaly Detection (EAD) Lateral Range Function 

generate or print a program data file (g/p)? 

By entering p, data can be printed from a program data file 
that was generated by the program. By entering g, a program 
data file can be generated for a set of user specified 
conditions. With either response, a sequence of additional 
queries is displayed. These queries require either an indicated 
response or a parameter value as the input. If the initial 
response is g, the sequence includes queries whose responses 
determine whether or not an auxiliary data file will be generated 
that can be used for future input of magnetic, processing or 
kinematic data. In particular, the first query in the sequence 
gives the option of using a combined electric, processing and 
kinematic data file. If the response indicates that it should 
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be used and the file is available, the parameter values that 
remain to be entered in order to generate a program data file are 
the following: the false alarm rate, the electric field noise, 
the maximum encounter lateral range and the lateral range step. 
The combined file should be used only if the effect of varying 
just one or more of these parameters is desired. If the response 
indicates that the file will not be used, queries concerning 
electric, kinematic and processing parameter values are 
displayed. 

After all of the program parameter values have been 
entered, a query is displayed giving the option of generating the 
combined file. Then a query giving the option of generating a 
program data file, a query giving the option of printing the 
encounter parameter values and a query giving the option of 
printing lateral range function values are displayed. The 
lateral range function values are the encounter detection 
probabilities indexed by lateral range. The parameters maximum 
lateral range and lateral range step determine the index lateral 
ranges of the encounters for which probabilities of detection are 
computed . 

The program generates electric signal values that correspond 
to points in time during an encounter. Following the lateral 
range function query, a query is displayed that gives the option 
of printing the electric signal values for an encounter. If the 
option is exercised, the option is repeated. When the option is 
not exercised, a query is displayed that gives the option of 
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generating or printing a new program data file. If this option 
is not exercised, the program ends. 

Some suggested guides for determining parameter values can 
be found in Section IV of this report. 
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Ill . Encounter Model Limitations 

In the two encounter models, the signal is the unfiltered 
signal that is generated by an electric dipole that moves 
relative to an electrometer. Depending on the input filtering, 
describing a submarine electric anomaly field as an electric 
dipole field should not be a significant limitation for an 
encounter slant range at the closest point of approach (CPA) that 
is greater than one hull length. The detection decision is based 
on samples from a single time interval (window) that is centered 
on the CPA. The length of the sampling interval and the sampling 
rate are parameter values that are inputs to the program. In 
terms of signal-to-noise ratio, there is an optimum sampling 
interval length (integration time) and sample rate. Although 
dipole signal energy is not symmetrically distributed about the 
CPA time, for a given sampling interval, the difference between 
the signal energy for the optimum interval location and the CPA 
centered location should not be significant in most cases. 

The encounter model electrometer noise samples are 
determined by a gaussian random process and they are values of 
independent identically distributed normal random variables. The 
standard deviation of these random variables is referred to as 
the electric noise and the variance is the noise in the sense of 
the signal-to-noise ratio. The degree of correspondence between 
this process and operational noise depends on the nature of the 
dominant operational noise sources and on the electrometer input 
filtering (noise whitening) . 
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In the encounter models, the intervals are adjacent but not 
overlapping and a detection statistic corresponds to each sample 
interval and its value is determined by the sample values. If 
the value of the detection statistic for an interval equals or 
exceeds a threshold value, a detection is indicated. The 
threshold value is determined by the false alarm probability 
which in turn is determined by the false alarm rate and the 
sampling interval length. 

A detector that used a moving sample interval that was 
generated by replacing the oldest sample by the newest one would 
correspond more closely to the detector in an operational 
detection system. In a model of the detector, since the sample 
windows overlap, the detection statistics would represent a 
sequence of random variables that were correlated over an 
interval equal to the width of the sample interval. Because of 
this dependence, it seems unlikely that the results that would be 
obtained with an encounter model based on an overlapping interval 
detector would differ significantly from those obtained with an 
encounter model based on a nonoverlapping interval detector. 

For encounter lengths of the order of a few nautical miles 
or less, the straight line encounter condition should not be a 
significant limitation. In particular, this should be the case 
for a fixed electrometer since, for a submarine target, vertical 
separation and course changes should be less likely to occur. 

Other models are available that can be used as the basis for 
computing an estimate of the probability that an electric anomaly 
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detection system will detect a submarine during an encounter. 

For example, one that is described in Appendix 6 could be used to 
determine the slant ranges of straight line encounters for which 
the detection probability is equal to a specified value. The 
parameter values that are required to do this are an average 
submarine electric dipole moment, a detection system capability 
factor and a noise factor. Values for these parameters can be 
determined from operational data. However, the values are 
specific to averages over a particular set of encounter 
conditions. An advantage of the two encounter models relative to 
this model is their adaptability to different electric, 
processing and kinematic conditions. 
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IV. Parameter Values 



The electric parameter value queries are generally 
explanatory with regard to the value that should be entered. 

This is also true of the kinematic parameter value queries. 
However, there is some ambiguity with respect the processing 
parameter value queries and the noise parameter value query. To 
reduce this ambiguity, a brief discussion of the common 
characteristics of the two encounter model processing parameters 
and the noise parameter is given below. This is followed by 
some guidelines for choosing these parameter values. 

In both encounter models, a decision is made at the end of 
each sampling interval. The decision is either noise energy was 
present during the interval or noise energy plus signal energy 
was present during the interval. The sample intervals are 
adjacent, equal width, nonoverlapping time intervals. The number 
of samples that are input in a sample interval is determined by 
the sampling rate and the interval length. 

The program default choice for the sampling rate is 2 • MAXF 
where MAXF is a parameter that is labeled the maximum electric 
signal frequency. This sample rate is the Nyquist rate for an 
ideal low pass filter. However, the signal in a sample interval 
that is computed by the program represents an unfiltered dipole 
signal. This is a reasonable approximation if the signal energy 
that is associated with signal components greater than MAXF is 
relatively small. As discussed below, the noise energy in a 
sample interval should be considered to be proportional to MAXF 
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in order to be consistent with the encounter models. Ideally, a 
default choice for MAXF would make the ratio of the signal energy 
to the noise energy a maximum for the sampling interval of an 
encounter. The program default choice for MAXF is 2 • MAXVM/MINRO 
where MAXVM is a user estimated maximum encounter relative speed 
converted from knots to meters per second and MINRO is a user 
estimated maximum slant range at CPA in meters in terms of a just 
detectable target. (A precise definition of a just detectable 
target can be made in terms of a specified detection probability, 
false alarm probability and target dipole moment.) The default 
choice for MAXF is consistent with the observation in Reference 1 
that if an optimum value for MAXF is determined for a minimum 
dipole moment target, then no significant increase in MAXF is 
required in order to maintain a required detection probability if 
the encounter lateral range is decreased even though the signal 
energy spectrum is shifted to higher frequencies. 

Because of the detection statistics and the gaussian noise 
model, if there were no penalty for decision delay, a sample 
interval length for a signal should be chosen equal to the signal 
duration, since this would make the detection probability for an 
encounter a maximum. In the program, a default choice for the 
sample interval length (integration time) is 2 • MAXRO/MINVM where 
MINVM is a user estimated maximum encounter relative speed 
converted from knots to meters per second and MAXRO is a user 
estimated maximum slant range at CPA in meters for a detectable 
encounter (in terms of a specific detection probability and false 
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alarm probability) with a user specified maximum dipole moment 
target. This choice might be considered a balance between 
minimizing decision delay and maximizing detection probability. 

In the encounter models, the sample intervals are located so that 
the CPA time is at the center of a sample interval and this 
sample interval is the only one that contains signal energy. 

This characteristic is consistent with the default choice for the 
sample interval. For a given sample interval length, although in 
general the interval location is not the optimum one in terms of 
signal energy, it should be approximately so in most cases. 

The program noise parameter is SIG. It represents the 
standard deviation o associated with the electric noise process 
of the two encounter models. In terms of the ideal low pass 
filter implied by the encounter models, its square should be 
equal to MAXF • (SIGO) s where SIGO is the electric noise process 
constant power spectral density. The program does not enforce 
this relation. Therefore, in using the program, the implied 
relation between the two input parameters: electric noise and 
maximum electric signal frequency should be kept in mind. If an 
average value of the peak-to-peak electric noise for an encounter 
can be estimated, for example from a electrometer trace, then the 
value for SIG should be chosen so that the estimate is 4 to 6 
times this value. 
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Appendix 1. The Detection Statistics 

In this appendix, yi, y 2 , • • • , y m are sequential values 

(voltages) representing the sample values in a sample interval. 
They are the input to an electrometer's detector. With these m 
sample values, the detector computes the value of a detection 
statistic. This value is represented by x and the detection 
decision corresponding to the sample interval is determined by 
the decision rule: If x > x*, then the input during the sample 
interval was noise plus signal, otherwise, the input was noise. 

For both encounter models, the detection probability and the 
false alarm probability are decreasing functions of x* and the 
relation is one-to-one in both cases. In the program, the false 
alarm probability Pf is used to determine a unique value of 
the threshold x*. This value is then used to determine a unique 
value of the detection probability p^. 

In the program, pf is found using the relation pf = RSt 
where R is the false alarm rate in false alarms per second 
and 6t is the sample interval length in seconds. This relation 
is based on the following argument: With no signal energy in a 
sample interval, yi = n lf y 2 = n 2 • ' ' ' < Yin = n m where ^ , n 2 , 

••• , n ro are noise values (voltages) input to an electrometer's 

detector. The nj_ are values of independent normal (gaussian) 
random variables, each with mean zero and standard deviation o. 
Because of this, in the encounter models, values of x for 
different sample intervals are the values of independent random 
variables that determine two outcomes: x > x* or x < x*. 
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Therefore, in terms of these outcomes, an encounter is a sequence 
of independent Bernoulli trials. If there is no signal present, 
since the noise is determined by a stationary process, pf is the 
same for each sample interval and, in terms of these outcomes, 
the sequence is a series of repeated independent Bernoulli trials 
and therefore the expected number of trials between false alarms 
is 1/Pf. Since the time between trials is St, the expected 
number of seconds between false alarm is St/pf or the expected 
number of false alarms per second R is equal to p f /St. 

The determination of x* depends on the encounter model 
statistic. For both encounter models, when there is a signal, 

Yl = n l + s l/ Y2 = n 2 + Y2> • • • / Ym = n m + s m where s 1# s 2 , 

• • • , s m are signal values (voltages) input to a electrometer's 
detector. The models imply that the signal values s^ = K- (E s ) ^ 
where K is a constant whose value is determined by the 
characteristics of the encounter electrometer and where the 
(E s )^ are dipole electric signal intensities. The models also 
imply that the noise values n^ are determined by a gaussian 
stochastic processes characterized by a standard deviation o 
and that n^ = K- (E N ) ± + n| where the (E N ) i are electric noise 
intensities and the n| are electrometer instrument noise values 
(voltages). In the program, the magnitude of K is 1. Since, 
for both models, p^ depends only on the ratio of signal energy 
to noise energy for a sample interval, this is a satisfactory 
choice for the program if instrument noise is assumed to be 
determined by a process independent of the electric noise process 
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and to be expressed in terms of an equivalent electric noise 
(E)i = (1/K) n- . 

For both encounter models, the signal (the average signal 

power) S = (1/m) -2 s^ and the noise (the expected value of the 

1 • . , 

average noise power) N = a so that the signal-to-noise ratio 

. 4 , 

is (1/m) -2 s i/a where the sum index l = 1, 2, • • • , m. 

The Crosscorrelation Detector Statistic: The statistic for the 

first encounter model is a crosscorrelation detector statistic 
that is defined by the sum 

x = 2 Yi-Sj. 

where the summation index i = 1, 2, • • • , m and the sum is over 

the values corresponding to a sample interval. For the first 
encounter model, the characteristics of both the noise and the 
signal are required in order to determine encounter detection 
probabilities. In particular, the signal values for an encounter 
are in the memory of the detector prior to the encounter. For 
the encounter conditions and a specified false alarm 
probability, the statistic is optimum in the sense that the 
encounter detection probability for this statistic is at least 
equal to that for any other statistic. Because of these 
considerations, encounter probabilities based on the 
crosscorrelation statistic can be considered to represent upper 
bounds on detection performance against dipole targets for 
electrometers of the type described by the models. 

For a sample interval without signal energy, x is the value 
of a normal random variable with a mean (i x = 0 and a variance 
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2 2 2 . . 
a x = o S where a is the standard deviation associated with 

the noise process and the sum index 1 = 1, 2, • • • , m and the 

sum is over the values corresponding to the sample interval. 

This implies that 

p f = 1 - P(x*/<?x) 

where P(z) is the standard normal cumulative distribution 
function. This relation is the basis for determining the 
threshold value x* . 

For the sample interval with signal energy, x is the value 

, 2 

of a normal random variable with a mean = s s i where the sum 

index i = 1, 2, • • • , m and the sum is over the values 

corresponding to the sample interval. This implies that 

Pd = 1 - p ( v * - dl ) 

where v* = x*/a^ and d = S s ±/o 2 = (1/m) -S/N. This relation 
is the basis for determining encounter detection probabilities 
for the first encounter model. The relation implies that for a 
specified false alarm probability pf the detection probability 
Pd is an increasing function of the signal to-noise ratio S/N. 
The Square Law Detector Statistic: The statistic for the second 

encounter model is a square law (energy) detector statistic that 
is defined by the sum 

x = z y[ 

where the sum index i = 1 , 2, • • • , m and the sum is over the 

values corresponding to the sample interval. For the second 
encounter model, only the characteristics of the noise are 
required to determine encounter detection probabilities. 
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is the 



For a sample interval without signal energy, x/a 2 
value of a chi-square random variable with m degrees of 
freedom. This implies that 

p f = 1 - P(x*/</ | m) 

where P(x*/o*|m) is the chi-square cumulative distribution 
function for a chi-square random variable with m degrees of 
freedom and where a is the standard deviation associated with 
the noise process. This relation is the basis for determining 
the threshold value x*. 

2 

For the sample interval with signal energy, x/a is the 

value of a noncentral chi-square random variable with m 

2 2 

degrees of freedom and noncentral parameter E s^/a where the 
sum index i = 1, 2, • • • , m and the sum is over the values 

corresponding to the sample interval. This implies that 

Pd = 1 “ P(x*/</|m, E s[/o 2 ) 

where P(x /a |m, E s^/o ) is the noncentral chi-square 

cumulative distribution function for a noncentral chi-square 

random variable with m degrees of freedom and noncentral 
2 2 

parameter E sj /a = (1/m) S/N. This relation is the basis for 
determining encounter detection probabilities for the second 
encounter model. The relation implies that for a specified false 
alarm probability pf, the detection probability p^ is an 
increasing function of the signal-to-noise ratio S/N. This is 
made more evident by the following relation: 

P(x*/a*|m, E s[/a ! ) = E{ (a^/ j ! ) • exp (-a) • P[x*/a * | (m + 2-j)]> 
where the sum index j = 0, 1, 2, • • • and the parameter 
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a = (1/2) E s l/o 2 = (m/2) • (S/N) where the sum index i 
• • • , m. Note that the relation > pf is satisfied 
P(x*/a J |m) < P[x*/a 2 |(m + 2j)] for j = 0, 1, 2, ••• 



= 1 , 2 , 
since 



17 



Appendix 2. Program Probability Calculations 

The program evaluates the cumulative and inverse cumulative 
distribution functions using approximations that are described in 
Reference 2. These approximations are listed below. 

The Standard Normal Cumulative Distribution Function 
Approximation : 

P(z) = 1 - s-t-(b]_ + t- (b 2 + t- (b 3 + t- (b 4 + t-b 5 )))) 
where s = ( 1/ (2 • ir ) * ) • exp (-z 2 / 2 ) and t = 1/(1 + b 0 z). And 
where 

b 0 = .2316419, b-L = .319381530, b 2 = -.356563782, 
b 4 = -1.821255978, b 5 = 1.330274429, 
and z > 0. For z < 0, P(z) = 1 - P(|z|). 

The Inverse Standard Normal Cumulative Distribution Function 
Approximation : 

Z (P) = t - (c 0 + t( Cl + tc 2 ) )/(l + t • (d x + t(d 2 + t-d 3 ))) 
where t = (ln(l/Q 2 ))5 and Q = 1 - P. And where 
Cq = 2.515517, c 1 = .802853, c 2 = .010328, 

d ± = 1.432788, d 2 = .189269, d 3 = .001308, 

and .5 < P < 1. For 0 < P < .5, z(P) = -z(l - P) . 

The Inverse Chi-Square Cumulative Distribution Function 
Approximation : 

v(p|m) = m • [ 1 - 2/ ( 9 • m) + z • (2/ (9 -m) ) § ] 3 

where P(z) = P(v|m). In the program, the inverse standard 

normal cumulative distribution function approximation is used to 
determine z. 
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The Noncentral Chi-Square Cumulative Distribution Function 
Approximation : 

P(w|m,E s i/o*) = P(z) 

where z= [2-w/(l + b)]$ -[2-a/(l+b)-l]| with 

2 2 2 2 2 2 
a = m + S Sj /a , b = (S s^/a )/(m + S s^/a ) and the sum index 

i = 1, 2, ••• , m. In the program, the standard normal 

cumulative distribution function approximation is used to 

determine P(z). 



19 



Appendix 3 . The Electric Signal 

The encounter models are defined by the following 
conditions: A submarine electric anomaly field is an electric 

dipole field. An electrometer electric signal value is the 
magnitude of the component of a dipole electric field in the 
direction of the electrometer axis. The basis for determining 
electric signal values in the program is an expression that can 
be developed as follows: In the right handed rectangular 
coordinate system that is shown in Figure 1, an electric dipole 
is at the origin, the xy-plane is the horizontal plane at a 
representative point in an encounter region, the positive y-axis 
is in the direction t of the electrometer's axis and the 
positive z-axis is positive upward. In this rectangular 
coordinate system, the unit vector with the orientation of the 
electrometer's axis can be expressed as e = j • cos $ + ksin $ 
where § is the electrometer axis depression angle. 

In a spherical coordinate system with the origin at the 
electric dipole and the polar axis in the direction of the dipole 
moment, = (c-p/r-*) • (2-r-cos © + 8- sin 0) is the electric 

field of the dipole at a point with spherical coordinates 
(r,r,6) . In this expression p is the magnitude of the dipole 
moment and c is a constant whose value is determined by the 
choice of units. The electric signal for an electrometer that is 
described by the encounter models is E s = e-E^ and it can be 
expressed in terms of the rectangular coordinate system as 
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z 




Figure 1. A unit vector in the direction of the electrometer axis 
and a unit vector in the direction of the dipole field 
of a dipole at the origin are shown at a point that is 
a distance r from the dipole. The unit vector h is 
in the direction of the horizontal component of the 
dipole moment. 
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follows: First, let r 0 be a unit vector with the orientation 

of the spherical coordinate system polar axis. Then the electric 

dipole field = (c-p/r 3 ) • (3 -r-cos e - r 0 ) since the unit 

1 

vector r= (i-x + j-y + k-z)/r where r = (x +y + z ) 5 , 

0 = [ (r 0 x r) x r]/sin 0 = [(r*r 0 )r - (r*r) • r 0 ]/sin 0, and the 
dot product r*ro = cos 0. The unit vector r 0 can be expressed 
in the rectangular coordinates by noting that r 0 = p/p and then 
expressing p in rectangular coordinates. To do this, let f2 
be the depression angle of p from the xy-plane (the horizontal 
plane) with ft positive downward and let a be the direction of 
p relative to the electrometer axis direction. Then the 
electric dipole moment p = p- (h-cos ft - k-sin ft) in terms of 
the unit vector k and the unit vector h = i-sin a + j • cos ft 
which has the direction of the horizontal component of p 
relative to the y-axis. Therefore, the unit vector in the dipole 
direction r 0 = i- (cos ft- sin a) + j • (cos ft- cos a) - ksin ft and 
E s = (c • p/r 3 ) • [ ( 3/r) • cos 0-(y-cos $ - z-sin §) 

(cos $-cos ft- cos a + sin f>-sin ft)] 
where cos © = (l/r)-(x-cos ft-sin a + y-cos ft-cos a - z-sin ft) 
since cos 0 = r*r 0 . As can be seen from this expression, for a 
constant dipole moment magnitude and direction and a constant 
electrometer direction, the electrometer signal is a function of 
the rectangular coordinates of the location of the electrometer 
relative to the dipole. In the encounter models, both of these 
conditions are satisfied. 
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Appendix 4 . The Anderson Formulation 

In the encounter models, the electric signal E s at a 
sample point in a straight line encounter can be represented in a 
form described by Anderson in Reference 4. For convenience, the 
Anderson formulation is used in the program to determine values 
for E s . It can be developed as follows: In Figure 2, the 

primed rectangular coordinate system is superimposed on the 
rectangular coordinate system of Figure 1 so that the origin is 
coincident with the origin of that system. An electrometer is in 
a straight line encounter with an electric dipole that is located 
at the origin of the combined system. The combined system moves 
with the electric dipole with the x'-axis oriented so that it is 
parallel to and in the direction of the electrometer's track 
relative to the dipole and the z'-axis oriented so that it is 
directed toward and passes through the CPA on that track. In the 
primed coordinate system, let 1, m, and n be the direction 
cosines of the dipole moment p and 1 1; and n! be the 

direction cosines of the unit vector e. Then, the unit vector 
r 0 = i' -1 + j' n + k' -m. And, for points on the relative track, 
x' = s', y' =0 and z' = R where R is the slant range of the 
dipole at CPA and s' is the electrometer's algebraic distance 
from CPA on the relative track. (It is negative for points before 
CPA and positive for points after CPA.) This implies that the 
unit vector r = i' • (s'/r) + k' • (R/r) . From Appendix 3, the 
field is = (c • p/r^ ) • (3 • r • cos 0 - rg) . With r and rg 

expressed in terms of the primed unit vectors and 
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Figure 2 . The dipole is at the origin in the unpriined and primed 
coordinate systems. In the primed coordinate system, 
the CPA is at (0, 0, R) and, for a time t during 
the encounter, the electrometer is at [s'(t), 0, R] . 
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